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Abstract — We focus on a particular form of network coding, 
reverse carpooling, in a wireless network where the potentially 
coded transmitted messages are to be decoded immediately upon 
reception. The network is fixed and known, and the system 
performance is measured in terms of the number of wireless 
broadcasts required to meet multiple unicast demands. Motivated 
by the structure of the coding scheme, we formulate the problem 
as a linear program by introducing a flow variable for each 
triple of connected nodes. This allows us to have a formulation 
polynomial in the number of nodes. Using dual decomposition 
and projected subgradient method, we present a decentralized 
algorithm to obtain optimal routing schemes in presence of coding 
opportunities. We show that the primal sub-problem can be 
expressed as a shortest path problem on an edge-graph, and the 
proposed algorithm requires each node to exchange information 
only with its neighbors. 

Keywords: Network coding, Distributed optimization. 



I. Introduction 

The network coding literature has provided a new set of 
tools for designing wireless networks. The early works in 
this literature were mainly focused on the multicast problem, 
where all sinks in the network require all of the information 
transmitted at the source node. Linear network coding is 
sufficient to achieve the capacity of a network for a multicast 
session (T). Another interesting problem is the multiple unicast 
problem, in which each source delivers its information to 
the corresponding destination. In this scenario, the coding 
schemes could be more complicated. In fact, linear network 
coding is not sufficient to achieve the capacity of a network 
with multiple unicast sessions [2|. Nevertheless, simple linear 
network coding techniques can boost the performance of 
the network significantly. An example of a simple coding 
opportunity across multiple sessions is reverse carpooling (see, 
for example [3]). Consider the three-node wireless network 
configuration depicted in Figure Q] where nodes A and B can 
only communicate via the relay node (R). The side nodes 
attempt to exchange messages Pa and Pb- Routing schemes 
require two transmissions at the relay node to exchange the 
messages. However, by allowing network coding we can save 
one transmission as follows. The relay node first receives 
the packets Pa and Pb from A and B, respectively. Then 
it forms the bitwise XOR of the packets and broadcasts the 
coded packet. Consequently, node A (B) can decode packet Pb 
(Pa) by XORing the received coded packet with the packet 
Pa (Pb) that it transmitted. This technique can be used to 
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Fig. 1. (a) Relay node requires two transmissions without network coding, 
(b) Relay node requires a single broadcast with network coding 



decrease the number of wireless transmissions (energy) when 
two unicast flows share a common path in opposite directions. 
Reverse carpooling is not the only form of opportunistic 
coding strategies. Katti et. al. JU present COPE, a practical 
scheme that exploits the overheard packets to introduce new 
coding opportunities. Experimental results in [4| demonstrate 
that COPE significantly increases the network throughput. 
Note that the broadcast nature of the wireless medium is the 
essential ingredient for such coding opportunities. 

In this work, we study the multiple unicast problem over 
a wireless network where only routing and reverse carpooling 
are allowed. The goal is to minimize the total cost of trans- 
missions (energy expenditure) to support the unicast sessions. 

Given the decentralized nature of the information, we are 
interested in distributed algorithms that can implement the 
optimal policies using coordination only among neighboring 
nodes. Note that without coding, this problem is equivalent 
to finding the shortest path from each source to the corre- 
sponding destination, which can be solved efficiently in a 
distributed manner by storing a routing table at each node 
and updating each table via local communication. Here, the 
routing decisions for each unicast session is independent of 
those for other sessions. In contrast, when reverse carpooling is 
allowed, the best path choice depends on the existence and rate 
of other unicast sessions. For example, consider the two unit- 
rate unicast sessions illustrated in Figure|2ja). In this example, 
deviation from the shortest path increases coding opportunities 
and hence decreases the total cost. Now, suppose the rate of 
session 2 is increased while the rate of session 1 is fixed. In this 
case, only session 1 deviates from its shortest path to reverse 
carpool with session 2. In this paper, using dual decomposition 
methods [8 |, we design a pricing mechanism at each node that 
allows us to obtain optimal routing decisions. Moreover, we 
show that optimal prices (dual variables) and routing decisions 
(primal variables) can be computed efficiently by exchanging 
information solely with immediate neighbors. 

Related literature addresses the opportunistic coding prob- 
lem from different points of view. Marden and Effros [5 | use 
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Fig. 2. (a) Both unit-rate unicast sessions deviate from the shortest path 
(direct line) to maximize coding opportunities, (b) The optimal routing 
decisions change when rate of session 2 is increased. 



a game theoretic approach which has become popular in the 
recent literature for devising distributed control policies for 
networks. This entails viewing each source destination pair as 
a selfish agent and designing a utility function for each one. It 
is well-known that such a multi-agent model creates incentive 
issues and leads to loss in performance. The authors of 121 
characterize suboptimality of the solutions reached by acting 
selfishly compared to the globally optimal solution. In contrast, 
we obtain the globally optimal solution using only limited 
coordination among the nodes. Also related to our work is 
Sengupta et al. H, which take a centralized optimization 
approach and provide various linear optimization formulations 
for a network with restricted forms of network coding. The 
formulations presented in (6J are based on all paths between 
each source destination pair, which lead to exponential-sized 
problems. In another related work, Reddy et al. Q adopt a 
path-based formulation together with population game dynam- 
ics to achieve a close approximation of the optimal policies. 
Our problem formulation, in contrast, is polynomial in the 
number of nodes of the network, and hence we are able to 
obtain the globally optimal solution in polynomial time. 

II. System Model 

Consider a wireless network with n nodes modeled as an 
undirected graph G = (Af,£), where N = {1, ...,n}, and 
{i,j} £ £ if and only if i and j can directly communicate to 
each other. We refer to such nodes as neighbors. We assume 
that the broadcast constraint is present, i.e, each transmission 
by a particular node is received by all of its neighbors. We also 
assume that the interference at the receivers is avoided by a 
Medium Access Control (MAC) protocol such as CSMA. 

We study a multiple unicast scenario where there is a set 
T of unicast flows. Each unicast flow i e T corresponds to 
a source-destination pair (st,dt), and a demand R t which is 
the rate at which packets from source s t should be delivered 
to destination d t reliably. Throughout this work we assume 
the set of unicast connections are feasible, i.e., there exists at 
least one path connecting each source to its destination, and 
no capacity constraint exists. 



For each node i, the cost of broadcasting a unit-rate flow 
(one packet per unit time) is denoted by c L . This could reflect 
the energy expenditure or the level of contention at a particular 
node. Throughout this work, we assume that only routing and 
reverse carpooling coding operations (as described in Section 
are allowed. The goal is to minimize the cost required 
to support all unicast connections in steady state. Next, we 
present a linear programming formulation for this problem. 

III. Problem Formulation 

In this section, we discuss how to formulate the prob- 
lem of coding-aware routing for multiple unicast sessions 
as a linear optimization problem. In general, there are two 
approaches to formulate a multiple-unicast problem which 
is also referred to as multi-commodity flow problem in the 
network optimization literature: The path-based approach and 
the link-based approach. The path-based approach considers 
all possible paths connecting each source-destination pair. This 
approach in general leads to a formulation of exponential 
size. Moreover, this approach is not suitable for distributing 
the decision making among individual wireless nodes. On 
the other hand, the link-based approach is more appealing 
from decentralized optimization point of view as well as the 
problem size. However, when reverse carpooling is allowed, 
there is no simple way of writing the transmission cost of a 
particular node in terms of the flow variables on its incoming 
and outgoing links. That is because in the case of reverse 
carpooling there is no one-to-one correspondence between 
packet transfers and physical transmissions. 

In order to have a finer description of the flows, we define 
the decision variables to describe the flow of packets over 
two successive hops. Let i be a node with neighbors v and w. 
Denote by x l (v,i, w) the flow of packets of the session t 6 T, 
transferred from v to w via node i. When reverse carpooling 
is allowed, node i can transfer J^teT xt ( v i h w ) packets from 
v to w, and J^teT xt ( w > h v ) packets from mtow with 

teT tGT 

transmissions. Since all reverse carpooling opportunities are 
along such triples of nodes, the total number of transmissions 
of a particular relay node i of graph G is 

where T G = j(i,j',fc) : £ £, {j, k} e£,i^ fcj. 

Note that every node can act as a relay, source or destination 
node. In order to write the constraints in a compact way, we 
define an expanded graph in which we introduce an artificial 
source and destination node for each unicast session. 

Definition 1. Given the network connectivity graph G = 
(TV, £), and a set of unicast connections described as 
{(s t ,d t , Rt) ■ t € T}, the expanded graph is an undirected 
graph G' = (Af', £') where 
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Fig. 3. (a) A graph G = (AT, €) with two source-destination pairs (b) 
Expanded graph G' = (Af',£') with artificial nodes added for each source 
and destination. 



£' = £ U (|J {{s' t ,s t },{d' t ,d t }]\ 

^ teT ' 

and the set of unicast connections is given by {(s' t ,d' t ,Rt) ■ 
teT}. 

Figure [3] illustrates a four node graph with two unicast con- 
nections (a) and its corresponding expanded graph (b). Every 
node of the expanded graph is either a source, a destination or 
a relay node. Moreover, all source and destination nodes of the 
expanded graph have degree one. This allows us to specify the 
conservation of flow constraints as follows: for every ( £ T 
and every ordered pair with {i,j} 6 £', write 



w.(i,j,w)er c 



E 



x t (v,i,j) = a t (ij), (2) 



Fig. 4. Conservation of flow constraints. The arrows demonstrate the direction 
of the flow of packets. 



the destinations of the expanded graph, one can write the total 
transmission cost as 

^ c i z l - y^ c dt R t , 
ieM teT 

where z% is the number of physical transmissions of node i 
on the expanded graph (cf. (dJ). Since, an additive constant 
does not change the optimal solution, we neglect the cost of 
retransmission at the destinations. 

We can now express the problem of minimizing the total 
transmission cost for a multiple unicast problem with reverse 
carpooling as the following linear program: 

minimize ^ E CiV ™ ^ 

ieAf v<w:(v,i,w)eT G / 

subject to 

^2x*(v,i,w) < y l vw , for all (v,i,w) E T G ,, 



where 



r G , = {(», j, k) : {*, j} e £', {j, k} e£',i? k} 



(3) 



is the set of all triples of connected nodes on the expanded 
graph G', and 



R t , if i = s' 
a t (i,j) = { -Rt, if j = d' t ; 

0, otherwise. 

The above constraints confirm that for any unicast session 
t, and all {i,j} G £, in each unit time, the difference between 
the number of packets that travel through node i to get to 
node j and the number of packets that travel from node i 
through node j to get to some other node equals either the 
number of packets that originate at node i if i is the source 
of unicast t or the number of packets demanded by node j if 
j is the destination of unicast t or zero otherwise (see Figure 
@J. For simplicity of notation, denote by the set of all 
non-negative vectors x f satisfying the above conservation of 
flow constraints. 

In order to correctly map the transmission cost on the 
expanded graph to that of the original graph note that each 
physical transmission of a particular node of the graph G maps 
to one transmission of the same node on the expanded graph 
G' . However, transmissions of the packets in the expanded 
graph that result in delivery of those packets to a destination 
do not correspond to any physical transmission in the original 
graph. By subtracting the cost of the extra transmissions to 



teT 



Vvw = Vlv, for a11 (v,i,w) g re, 
a;'eA (t) , for all teT, 

where Tc is defined in (0), and Aw is the set of non-negative 
vectors x* satisfying the conservation of flow constraints 
(ffji. The preceding formulation of the problem is a linear 
optimization problem of polynomial size which can be solved 
efficiently in a centralized manner. However, we are interested 
in a distributed algorithm to solve this problem. In the fol- 
lowing section, we use a dual decomposition together with 
subgradient method to solve this problem in a decentralized 
manner. 

IV. Dual Decomposition 

We consider the following Lagrangian function that penal- 
izes the first set of constraints of the problem in (|4j. For each 
triple (v,i,w) £ Tg>, the Lagrange multiplier corresponding 
to the constraint for (v,i,w) is denoted by p(v,i,w). 



L{x,y,p) = Y E ! ' // - 

ieJV v<w.(v,i,w)EV c 



teT 



Vwv 



The dual function is thus given by 
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q(p) = min 



E E 



E E 

isW i)<«i,(u,i,ii))6r G i 

subject to 



p(v, i,w)x (v, i,w) (5) 
Vlw ( c i ~ p( v > h w) - p(w, i, v)) 



Vvw Vwv> 



for all (v, i, w) G Tc, 
for all t G T. 



Since y l v can be positive or negative, we can drive the cost 
to — oo when Ci — p(v,i,w) — p(w,i,v) is non-zero for any 
triple (v,i,w). Hence, the dual function can be simplified as 

!q*(p), if (cj -p(v,i,w) -p(w,i,v)) = 0, 
for all (v, i, w) G T G r 7 
— oo, otherwise, 

where q*(p) is the optimal value of the primal sub-problem 
given by 

g*(p)=min p(v,i,w)x t (v,i,w) (6) 

subject to 



x L G 



for all t G T, 



(7) 



The t/Mfl/ problem solves the maximization 

maximize q* (p) 
subject to 

p(v, i, w) + p(w, i, v) = Ci, for all (v, i, w) G Tc, 

p(v, i, w) > 0, for all (y, i, w) G Tq 1 - 

Note that for every feasible vector p, q* (p) provides a lower 
bound on the optimal cost of the problem in (0J. Since there is 
no duality gap for feasible linear optimization problems (see 
||9l), optimizing the dual problem achieves the optimal value 
of the original problem (0). 

We next show how to solve the dual problem given an opti- 
mal solution to the primal sub-problem and how to decompose 
and solve the primal sub-problem efficiently. 

A. Dual Problem 

We employ the projected subgradient optimization method 
(cf. Section 6.3.1 of [9|) to solve the dual problem given by 
<j7j. We start from a dual feasible solution p[0}. Given the n th 
iteration p[n] for any n > 0, we compute a subgradient of 
q(p) at p[n], which is given by 



where y[n]) is an optimal solution of the problem 

in (|5]l for the dual vector p — p[n]. Note that for a dual 
feasible solution, the values of the y l vw do not affect the cost 
or feasibility of problem @. Hence, they can be selected 
arbitrarily, and we choose them to be zero. Therefore, x[n] 
can be computed by solving the primal sub-problem The 
(n + l) st iterate of the subgradient method is given by 

p[n + l] = (p[n]+ a n ^ sfe*[n]) + , (8) 

t<ET 



x(w— *v) 




Fig. 5. An iteration of the projected subgradient method. x(v — ¥ w) = 
a n 5Z tgT i, w)[n] is proportional to the total flow relayed from v to 
to through i. Observe that when x(v —>«>)> x(w — > v) the price of 
forwarding packets from w to v, p(w, i, v), decreases to attract more flow in 
this direction. 



where a n is the stepsize, and (■)+ denotes Euclidean projec- 
tion on the set of feasible solutions of the dual problem, i.e., 
z + is the optimal solution of the following problem 

min 2_. (z(v,i,w) — p(v,i,w)) 2 (9) 

2 



-{z(w,i,v) -p(w,i,v)Y 



subject to 



p{v, i, w) + p(w, i, v) — Ci, for all (v, i, w) G Tc, 
p(v,i,w) > 0, for all (v,i,w) G r<y. 

We can decompose the above projection problem into 
smaller sub-problems of projection on a two-dimensional 
simplex, which can be expressed in closed-form (see Figure [5] 
for an illustration of the process). After a few straightforward 
steps, we can verify that for all triples (v,i,w) G Tq' 



p(y, i, w)[n + 1] = p(v,i,w)[n] 

(x* (v,i,w) [n] — i* (w, i,v)[i 



a" 

"T 



E 



[0,c 4 ] 



(10) 




where [z]^ a b j denotes projection of z into interval [a, b], i.e., 

if z < a 
if a < z < b 
if z > b. 

Under proper stepsize selection rules, such as the diminish- 
ing stepsize selection rule a n = —, the iterates given in (fTOb 
converge to the dual optimal solution p* as n — > oo [9|. Even 
though the optimal solution of the primal sub-problem © at 
p* gives the optimal solution of the main problem ©, there 
is no guarantee that the corresponding primal solutions x[n] 
approach the optimal solution of the main problem © as p[n] 
approach p* . However, there are several methods suggested 
by Sherali and Choi IfTUl to recover close to optimal primal 
solutions from the sequence of primal iterates {x[n]}. For 
instance, consider the sequence {a;[n]}, where 



fe=l 



x[k] 



(11) 



For proper stepsize selection rules, any limit point of the 
sequence gives an optimal solution to the main problem 

© by UTOj. 
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B. Primal Sub-Problem 

We proceed to the primal sub-problem described in ©■ 
Since the cost function is separable and there is no constraint 
binding the flow variables i, w) for different t G T, we 
can decompose this problem into separate sub-problems. Write 

q*(p) = X^teT S*(p)> wnere f° r eacn t ET 

<7 t (p)=min p(v,i,w)x t (v,i,w) (12) 

(u,i,u;)£r G / 

subject to 

^ x t (i,j,w) - ^2 x t (v,i,j)=a t (i,j), 

w:(i,j,w)er G , v:(v,i,j)eF a i 

for all s.t. {i,j} G 
i, w) > 0, for all (v, i, w) G Tc. 

Each such problem is a minimum cost flow problem with 
no capacity constraint. We show that this problem can be 
expressed as a shortest path problem on an edge-graph defined 
below. 

Definition 2. The edge-graph H = (V, A) for the undirected 
graph G = {JV,£) is defined as a weighted, directed graph 
where 

V = Ui,j):i,jeAr, wd{i,j}eS'\, 
A = {((*,j),(j,fc)) : (i,j,k)eT c 

In other words, there is a directed link from node G V 
to (fc, Z) G V if and only if j = k, and i ^ I, 

An example of an expanded graph with one source- 
destination pair and the corresponding edge-graph is illustrated 
in Figure [6] Next, we establish the equivalence relation. 

Proposition 1. For each t G T, the problem in ( T/21 ) is 
equivalent to the shortest path problem between nodes (sj, St) 
and (d' t , d%) of the edge-graph H defined in Definition^ when 
the cost of each link (j, fc)J G A is given by p(i,j, k). 

Proof: Consider the problem of finding the shortest path 
between nodes p and q of a directed graph (V, -4), when the 
cost of each link G A is denoted by Wij. A linear 

programming formulation of this problem is 



minimize WijXij (13) 

subject to 

Xij — xji — Oi for all i G V, 

> 0, for all G A, 

where 

1 , if i = p ; 
— 1, if i = q; 
0, otherwise. 

In order to see the equivalence of this problem to the 
problem in ( fT2l . it is sufficient to use the construction of 
the edge-graph H given by Definition [2] with p = (s' t ,St) 




(a) 



(b) 



Fig. 6. (a) Expanded graph with source node s', and destination node d' . 
(b) Corresponding edge-graph generated as described in with (s',1) as 
the source and (3, d') as the destination. The solid links can participate in a 
path from the source to the destination. 



and q = (d' t ,d t ). Note that there is a one to one cor- 
respondence between triples (v,i,w) G Tq>, and directed 
links [(v,i), (i,w)) G A, and the weight of each such triple 
coincides with the cost of the related link. We also need to 
scale the variables by R t to obtain the exact same formulation 
as in (TPH l. but this does not affect the equivalence. ■ 
There are various decentralized and asynchronous imple- 
mentations of the shortest path problem available in the 
literature such as the well-known Bellman-Ford and Dijkstra's 
algorithms, (see [8 |). However, distributed computation of the 
shortest path on the edge-graph does not immediately yield 
a decentralized algorithm on the original network of wireless 
nodes. We discuss distributed implementation of the presented 
algorithms on the original graph in the following section. 

V. Decentralized Implementation 

We provided a linear programming formulation of a mini- 
mum cost coding-aware routing problem for multiple unicast 
sessions. We presented an iterative algorithm based on the 
projected subgradient method to solve this problem. Each 
iteration of the algorithm is computed as in ( TTOb given an 
optimal solution of the primal sub-problem in ©. We now 
discuss how to implement this iterative algorithm when nodes 
can only exchange information with their neighbors. 

The dual (price) variables p(v,i,w) as well as all primal 
(flow) variables of the form x*{v, i,w) are stored at node 
i. Given an optimal solution to the primal sub-problem, the 
subgradient iteration given in ( TTOb can be readily computed at 
node i without any further coordination with other nodes. 

Now consider the primal sub-problem ©. As we discussed 
earlier, the primal sub-problem is equivalent to a shortest 
path problem on the edge-graph described in Definition |2] 
for each source-destination pair. The shortest path problem on 
the edge-graph can be solved by the Bellman-Ford algorithm 
in a decentralized and asynchronous manner (cf. J8)). In this 
algorithm, every node (i,w) G V updates its label based on 
the information obtained from its predecessor (v,i) G V, as 
well as the cost of the link between them given by p(v, i, w). 

For every node of the edge-graph H, assign its 

processor to node i of the expanded graph G' , i.e., the distance 
information kept at node of the edge-graph is stored at 
node i. Hence, information exchange between nodes (i, w) and 
(v, i) in the edge-graph is equivalent to information exchange 



Fig. 7. Network connectivity graph. 

between nodes i and v in the expanded graph. This is possible 
because nodes i and v of the expanded graph are connected 
by definition of V in Definition [2] Therefore, the primal 
sub-problem §6^ can be solved on the expanded graph by 
exchanging information only with immediate neighbors. This 
immediately yields a distributed algorithm on the original 
graph G by performing the computations assigned to each 
source (destination) node of graph G' at the corresponding 
source (destination) node of the original graph G. 

We can interpret the operation of the algorithm as follows. 
At each round, each node i declares a price p(v, i,w) for 
forwarding unit flow from v to w, and a different price 
p(w,i,v) for the opposite direction. Given these prices, the 
nodes obtain the shortest path from each source to the cor- 
responding destination. This gives the updated flow variables, 
i.e., how much flow i needs to relay from vtow and how much 
is relayed from w to v. Finally, node i updates the prices so 
that the price of forwarding packets in the direction with less 
flow decreases (see Figure |5}. This provides an incentive for 
other flows to exploit more coding opportunities. 

VI. Simulation Results 

We demonstrate the convergence of the presented iterative 
algorithm for solving the minimum cost coding-aware routing 
problem (|4j. The connectivity graph G is generated by ran- 
domly placing nodes in a square according to a Poisson point 
process of unit rate. Two nodes are assumed to be connected 
if and only if their distance is less than one. Figure Q illus- 
trates a sample graph with the following source-destination 
pairs: (s u di) = (20, 13); (s 2 , <h) = (26, 7); (s 3 , d 3 ) = 
(15, 23); (s4, di) = (7, 22). We assume that each pair has unit 
rate, and the cost of each transmission is also one. 

We start the subgradient iterations in (TlOb with the initial 
dual variables given by p(v, i,w) = ^. We take a diminishing 
stepsize selection rule a™ = ~ and recover the primal feasible 
solutions x[n] as in ( fTTT i. Figure [8] shows the total system cost 
evaluated at x[n] for each iteration of the subgradient method. 
The optimal value of the primal sub-problem (O at any dual 
feasible solution provides a lower bound on the optimal cost 
which is plotted in Figure [S] The total transmission cost for a 
plain routing scheme is also plotted as a reference. Note that 
after a few iterations of the subgradient method the transmis- 
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Fig. 8. Convergence of the iterations of the subgradient method. 

sion cost of the system with reverse carpooling significantly 
improves compared to that of plain routing schemes. 

VII. Conclusions 

We studied the minimum cost multiple unicast problem over 
a wireless network when coding opportunities of the form 
of reverse carpooling exist. We formulated the problem as a 
linear optimization problem of polynomial size by considering 
flow variables across triples of connected nodes. By employing 
dual decomposition and the projected subgradient method, we 
provided an iterative algorithm to obtain the optimal flow 
control strategies. We also showed that the presented algorithm 
can be implemented in a distributed manner so that each node 
only requires to exchange information with its neighbors. 

The problem formulation in this work as well as the ones 
in (6) and Q are based on maximizing coding opportunities 
to decrease the total transmission cost. However, this can 
cause a severe congestion in a certain area of the network, 
while there are other unused paths that can provide similar 
coding opportunities. In such scenarios, it is more desirable 
to distribute the flow evenly across multiple paths to avoid 
congestion. We shall study this problem in the future works. 
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